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^ : Abstract 

We review the occurrence of electric-field domains in doped superlattices 
' within a discrete drift model. A complete analysis of the construction and 

stability of stationary field profiles having two domains is carried out. As 



> 

in 



a consequence, we can provide a simple analytical estimation for the doping 
density above which stable stable domains occur. This bound may be useful 
Tlj" I for the design of superlattices exhibiting self-sustained current oscillations. 

I Furthermore we explain why stable domains occur in superlattices in contrast 

O ' to the usual Gunn diode. 

^ : I. INTRODUCTION 

a 

^ I Electrical transport in semiconductor superlattices (SL) has attracted much properties 

pH ■ related to the artificial band structure. One of these features is the occurrence of stationary 
P ■ electric field domains which have already been reported in Ref . . Due to advanced growing 
facilities and experimental techniques the complicated structure of the current- volt age char- 
acteristics which exhibits several branches, roughly equal to the number of quantum wells, 
^ I could be resolved during the last decade 0-0]. In these experiments it was demonstrated 
■ that the different branches are connected to the formation of two domains with different 
values of the electric field inside the sample. Depending on the conditions, stable stationary 
field domains and traveling domain boundaries may occur. In the latter case, the dynamics 
of the electric field domains gives rise to time-dependent oscillations of the current. P|-[TT[] 

Early approaches towards a modelling of these phenomena have been performed within 
a continuum model including strong diffusion |T^ and with the help of an equivalent circuit 



|13| . During the last years it has been shown theoretically that the observed phenomena can 
be fairly well reproduced by models which essentially combine the discrete Poisson-equation 
and rate equations for the carrier densities in the different quantum wells |I4|-[T6[]. Also 



the time-dependent current oscillations could be recovered in these models [|T^,|l^,|l3] • ^ 
prediction of spatio-temporal chaos in resonant-tunneling superlattices under dc+ac voltage 
bias has been made on the basis of the discrete drift model The influence of growth- 
related imperfections on the SL behavior has been studied in [^]. Some of these phenomena 
may also be described by discrete models with Monte Carlo dynamics incorporating single- 
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electron tunneling effects (wich are important for slim superlattices and give rise to additional 
oscillations of the current ||2l|j2^). 

In this paper we want to explain how these complicated phenomena are generated by 
such models. This provides a deeper insight into the basic mechanisms and helps to classify 
the results of various experiments and computer simulations. In particular we want to give 
an answer to the following questions: How is it possible to understand the appearance of 
the complicated structure of the current- voltage characteristic? What are the conditions for 
stability and oscillations and how can they be understood? What is the main difference to 
the Gunn diode, where hardly any stable domain states are observed? 

The paper is organized as follows: The model we use is described in the second section. 
The third section shows how the complex stationary current-voltage characteristic changes 
as the doping increases. In the fourth section we investigate the stability of the stationary 
states and prove an explicit criterion for the occurrence of stable domain states. The last 
section contains our conclusions and the Appendix is devoted to a proof that no self-sustained 
oscillations appear for a small product of doping times the number of SL periods. When the 
doping density lies in the range between these two stability boundaries, the model exhibits 
self-sustained oscillations. Results concerning self-sustained oscillations of the current are 



presented in a companion paper |]23|, where a direct comparison with experimental data is 
made. 



II. THE MODEL 

We consider a semiconductor superlattice where the lateral extension of wells and barriers 
is much larger than the total length of the SL, so that single-electron tunneling effects (see e.g. 
Ref. ||2T|) are negligible. The quantum wells (QW) are weakly coupled and the scattering 
times are much shorter than the tunneling time between adjacent QWs. Thus it makes 
sense to consider the electrons to be localized within the QWs and in local equilibrium at 
the lattice temperature. The current is mainly determined by the resonances between the 
different energy levels in the QWs, which we denote by Ci, z = 1, 2, . . ., in order of increasing 
energy counted from the bottom of the conduction band. For the biases of interest here, 
there are three important resonances ClCl, C1C2 and C1C3. If the intersubband relaxation 
is also fast with respect to the tunneling, in practice only the lowest subband is occupied. 
In this case it is a reasonable approximation to consider the QWs as entities characterized 
by average values of the electron density nj in the i-th QW (in units cm~^) and the electric 
field Ei between wells i and i -|- 1, with i = 1, . . . N. For sufficiently high electric fields, the 
electrons tunnel only in the forward direction and we can describe the dynamical behavior 
of hi by the rate equation 

^ = J (v{Ei^i) hi^i - v{Ei) hi^ , (1) 

where / is the superlattice period and v{E)/l is the average electron tunneling-rate for the 
local field. v{E) has a peak at certain values of the electric field connected to resonant 
tunneling ClCl, C1C2 and C1C3. In this paper we shall be concerned with phenomena 
occurring at fields higher than the first resonant peak C1C2, so that we shall omit the 
miniband peak ClCl in our tunneling rate v{E)] see the curve plotted in Fig. |l]by a dashed 
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line. Phenomena at lower fields have been studied by adding the ClCl peak to our curve, 



2J]. The electric fields must fulfil the Poisson equation averaged over one SL period 



Ei - E^-i = ^ (n,, - Nd) , (2) 

where e, g, and Nd are the average permittivity, the charge of the electron, and the donor 
density per SL-period (in units cm~^), respectively. By differentiating Eq. (0) with respect 
to time we obtain Ampere's law for the total current density J{t) [|17|1 : 

e ^ + lv{Ei) (qND + e{% - = J, (3) 

at I ^ ' 

where i = 1,. . . ,N. Typically the experiments are performed with a constant dc voltage 
bias $ yielding 

N _ 

lY.E^=~^. (4) 

i=l 

Notice that there are 2N+2 unknowns: 

Eq, -El, . . . , Em-, hi, ... , un, J 

and 2N+1 equations, so that we need to specify one boundary condition for £"0 plus an 
appropriate initial profile Ei{0). The boundary condition for Eq (the average electric field 
before the SL) can be fixed by specifying the electron density at the first site, hi, according 
to d^). In typical experiments the region before the SL has an excess of electrons due to a 
stronger n-doping there than in the SL [|8|-[TT|. Thus it is plausible assuming that there is 
an excess number of electrons at the first SL period measured by a dimensionless parameter 
-1: 

hi = {l + c)Nd ■ (5) 

c has to be quite small because it is known that a steady uniform-electric-field profile is 
observed at low laser illumination in undoped SL pJT^Jl^. This observation allows us to infer 



the electron tunneling-rate v{E) directly from measured current- voltage data Another 
possibility is to derive v{E) from simple one-dimensional quantum- mechanical calculations 
of resonant tunneling, as was done by Prengel et al IQ. They used a more complicated 
discrete model with two electron densities corresponding to the populations of the two lower 
energy levels of each QW. Their model reduces to a form of ours when the large separation 
between the time scales of phonon scattering, resonant tunneling and dielectric relaxation 
is taken into account. 

For the calculations that follow, it is convenient to render the equations (0)-(^ dimen- 
sionless by adopting as the units of electric field and tunneling-rate the values at the C1C2 
peak of the v{E) curve, v{E), Em and vm (about lO^V/cm and 10'^ cm/s, respectively, for 
the sample of Ref. |10|)- Like in Ref. we set: 

Ei qhi Jl 

Ei = , n 



Em ' ^ eEM ' qNdVm ' 

— , t = = — ^ , = — = — = (6) 

VM cEmI NEmI 
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We obtain the dimensionless equations: 



^ = m - (l + ^^^^) v{E^) (7) 

for i = 1, . . . iV, 

1 ^ 

<p=j^j:m, (8) 

and the boundary condition 

Eo{t) = E^{t) ~ cu. (9) 
The dimensionless parameter u, is defined by 

Nnq 



eE 



(10) 



M 



which yields about 0.1 for the SL used in the experiments [|TO|,|TT|). The constant voltage 
condition (H) determines the current to be 

1 ^ / E- - E' i\ 

m = i^E{^ + ^^^^^ (11) 

With the choice (|^), the dimensionless tunneling rate v{E) has a maximum at i^^ = 1 with 
f (1) = 1. Throughout this paper we use the function v{E) which is plotted in Fig. |l] by a 
dashed line. Besides having a maximum f (1) = 1, it has a minimum at Em ~ 1.667 with 
v{Em) = I'm ~ 0.323. Nevertheless nearly all of the features discussed in the following are 
independent of the exact shape of the function v{E). We only impose the restrictions that 
v{E) > for E > and the existence of a minimum at Em > 1. 

For any > 0, we shall assume that the initial electric field profile is strictly positive and 
that the electron density is non-negative: Ei{0) > 0, nj(0) = -E'j(O)— -E'i_i(0)+z/ > 0, Wi. This 
is reasonable unless (j) is close to zero (but then the ClCl peak in the tunneling rate curve 
should be restored) or the boundary conditions are unrealistic. From the equations and our 
assumption on the initial field profile, it follows that Ei > and rij = Ei — Ei_i + v > for all 
positive times. The model equations have interesting properties concerning the monotonic 
behavior of the electric fields with respect to the QW number. These are summarized in the 
following Lemmas: 

Lemma 1. // the fields of two adjacent QWs are identical, i.e., E^ = Ek_i holds for 
some k with 2 < k < N , there is at least one i from 1 < i < k with Ei = Ei_i and 
d{Ei — Ei^i)/dt 7^ 0. For c = there is the additional possibility that Eq = Ei = . . . = Ek 
holds. 



Lemma 2. // c > and the field distribution is monotone increasing at t = (Ei{0) > 
Ei-i{0), ^i), it will will keep this property for all later times t > 0. (For c < 0, the same 
holds for a decreasing field distribution). 
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Both lemmas can be easily proved. If Ek = E^-i holds Eq. ^ gives 

^{E, - E,_,) = Ehzl^^zl (12) 
at V 

Therefore diE^ — Ek-i)/dt = implies Ek^i = Ek-2- Then we may repeat the argument 
until we either find an i G {2, 3, . . . k} with Ei = Ei_i and d{Ei — Ei_i)/dt 7^ or we find 
El = Eq which violates the boundary condition for c 7^ 0. This proves lemma 1. Lemma 
2 is proved by contradiction. Assume that at t > the sequence {E^} is increasing but 
Ek = Ek_i and d{Ek — Ef,_i)/dt < for some k. Then Eq. (^) yields E^^i — Ek-2 < in 
contradiction to the assumed increasing behavior at t. 

In order to understand the properties of the current-voltage characteristics, we will use 
repeatedly these basic lemmas in what follows. 



III. STATIONARY STATES 



In this chapter we want to explain how the complex domain structures found experimen- 
tally and from computer simulations PJ^,^,^ are generated by this simple model. 

We denote the electric field profile and the current density of stationary states by E* and 
/*, respectively. An easy way to construct the stationary profiles is to fix /*, find out the 
corresponding electric field profile {E*}, i = 1, . . . , N, calculate their voltage as a function 
of/* 

1 ^ 

<^(n = T7E^r(n- (13) 

i=i 

The field profile must fulfil the equation 

Ell = E: + u(^1- =: f{E:, n (U) 

which was also used in Ref. [|1^. The boundary condition implies 

which has three solutions E^ for a known fixed value of the current on the interval (1+c) Vm < 
r <l + c. 

In order to understand the properties of the stationary profiles we will now investigate 
the behavior of the set {E*} as a function of E^. At first we restrict ourselves to c > and 
increasing field profiles. To construct {E*}, we have to invert the function f{E,I*) for a 
fixed value of /*. Its derivative is: 

dE v{Ey dE ^ ' 

With the restriction to increasing field profiles, we can always obtain E* for E^ > Em 
because dv/dE > and f{Em, I*) < Em holds. If E^ < Em, Eq. (|15D implies /* < 1 + c and 
we find that f{E, I*) is strictly monotone increasing for all E if 
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(17) 



yielding a condition which only depends on the v{E) relation but is not dependent on I. For 
our function v{E) this yields u < 0.195/ (1 + c). In this case the function f{E, I*) {I* fixed) 
is always invertible. Then we can find a unique field profile parametrized by the point E^. 
Since we have three possible solutions El of (|15D for each given /*/(l + c) e {vm, 1), there are 
three different voltages for each value of the current in this range. The function 0(/*) is 
thus three- valued, which means that by inverting it we obtain an N- or an Z-shaped current- 
voltage characteristic as shown in Fig. |l| for u = 0.05 and u = 0.15, respectively. Both types 
can be easily understood: When the doping density u is low, Eq. (p^ ) shows that E* ^ £"* 
holds. Thus, the field profile is nearly uniform and the current-voltage characteristics follows 
the v{E)-cmve as shown in Fig. |1| for u = 0.05. This is physically obvious as there are few 
charges present inside the sample. For larger values of u the values E* may strongly deviate 
from E^ with increasing i, ii E^ is not a fixed point of f{E,I*) (which is the case for c = 
1|1). Let us denote by E'-^^I) < E^'^\l) < E^^\l) the three solutions of v{E) = / for a 



given / G (fm, 1) which are the fixed points of Eq. (|T^. If is small and c > 0, E^ is on 
the first branch of v{E) and the values E* tend to E^^\l*). When is larger and El is 
located on the second branch of v{E), the sequence E* leaves the neig hborhood of E^^^I*) 
and then approaches E^^\l*) on the third branch of v{E) if u is large enough. This is shown 
in Fig. ^(a). In this case the voltage is basicaly determined by the fixed point E^^\l*). 
Since dE^^\l) / dl > 0, this branch of stationary solutions may exhibit a range of positive 
differential conductance leading to the Z-shape. This effect is more pronounced for longer 
superlattices with many wells N and also for larger values of c. 

For larger doping u the condition (|1^ is violated and the function f{E,I) may not 
be invertible for some current /. In this case there can be more than one possible -Ej+i 
following a given Ei. Then the current- voltage characteristic can no longer be unambiguously 
parametrized by the point El. In general f{E,I) has 3 different branches for a certain 
interval of /*, as shown in Fig. ^(b). Let us call branch a that having df /dE > for low 
E, branch (3 has df /dE < 0, and branch 7 again has df /dE > but for larger E. 

Let us explain how to construct different stationary field profiles for a given value of the 
current /*. We shall assume that the profiles are increasing, E*_^-^ > E*, Vz. First of all. El 
may be located on branch 7 of f{E, I*), and so will all successive fields E*. This profile will 
have the largest possible voltage for the same I*. Secondly, El may be located on branch 
f3, which implies that all successive E* of an increasing field profile have to be on branch 
7. The corresponding voltage is smaller than that of the previously described branch but 
larger than those stationary solution branches that we analyze next. 

If El is on branch a, we may have E* , i = 1, . . . , j — 1, (j = 2, . . . , N) on branch a, and 
E* either on branch (3 or on branch 7. We obtain a different branch of stationary solutions 
for each such possibility. Let us denote by (j, (3) or (j, 7) the solution branch having E* 
either on branch /? or 7 respectively, and E*, i = 1, . . . , j — 1 on branch a. In Fig. ^(b) a 
solution {j, 13) is shown by a dashed line and a solution (j, 7) by a full line. Clearly j = 1 
corresponds to the possibilities discussed above. Finally, we have one solution where all field 
values are on branch a which we denote by (A^ -t- 1,7). In order of increasing voltages, we 
have 



6 



0(7V+1,7)(^*) < (PiN,/3){I*) < 0(iV,7)(/*) < 
...<0(l,/3)(/*)<0(l,,)(r), 

corresponding to 2N + 1 different stationary solution branches with the same current J*. 
They can be observed in Fig. |l| for u = 1.0 and /* = 0.8. Notice that the branches (j + 1, 7) 
and coalesce at a current /* G (1, 1 + c) which is roughly independent of u (see Fig. 

H(c)). The branches and (j, 7) coalesce at a lower current Ic which decreases as u 

increases (see Fig. 0(d)). The current- voltage characteristic curve is thus connected as 
shown in Fig. |I| (i^ = 1.0). 

The field profile of the solution branch (15,7) is depicted by the crosses in Fig. |^. One 
can clearly identify two regions 1 < i < j and j < i < N where the electric field is 
roughly constant and close to a fixed point with v{E*) ^ I*. In between there is a transition 
layer, the domain boundary, consisting of only a few wells. These type of states we call 
domain states. A shift of the domain boundary by one well only changes the voltage as long 
as the transition layer does not extend to one of the contacts. As the stationary solutions 
resulting from a one-well shift are very similar, the domain branches in the current-voltage 
characteristics look alike, as can be seen in Fig. (z^ = 1.0). The slope of the different domain 
branches in Fig. may vary from branch to branch. The high field domain coresponding to 
far right branches has a much larger extension than the low field domain. Then the slope of 
these branches in Fig. |1] will be closer to the slope (conductivity) of the third branch of the 
curve v{E). Similarly the slope of far left branches in Fig. will be closer to the slope of 
the first branch of v{E). The larger the difference in slope between first and third branches 
of v{E) is, the larger the variation in the slope of the branches in Fig. [l| will be. 

Note that the domain states are not very sensitive to the exact type of boundary condi- 
tions if two conditions are fulfilled: 

1. The boundary conditions must allow for the existence of a roughly constant field 
distribution E* ^ E^-^^ and E* ^ E^^^ probably after a short contact layer of some 
wells. 

2. The domain boundary must be located sufficiently deep inside the sample, so that it 
does not collide with the contact layer. 

As u decreases, the solution branches become shorter and eventually disappear if Ic 
becomes larger than (1 + c), which happens if the inequality (0) holds. The stationary 
domain structures are seen for narrow current intervals about /* = 1 for intermediate doping 
as shown in Fig. |l|(z/ = 0.3). Another complex feature can be found here for larger voltages 
where extra wiggles appear. They occur if E^ crosses the value 1, yielding an additional 
maximum in I*. 

Except for the fundamental question of stability, we have now understood the morphology 
of the complicated current-voltage characteristic curve shown in Fig. |I], its changes with 
doping and its relation with the electric field profile. The very same features occur in the 



more complicated model of Prengel et al. p!^ , p!8| as shown numerically in Ref. p5 |. 

So far we have restricted ourselves to increasing field profiles. Therefore we could only 
find domain states where the high-field domain is located at the receiving contact and the 
domain boundary is an accumulation layer. Nevertheless, this is not the full story. For 
sufficiently large i> other solutions are also possible, even for c > 0. A typical such field 
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profile is depicted by circles in Fig. |^. The field starts on branch 7 of the function f{E, I*) 
and first increases with the QW index towards the third fixed point E^^\l*). At a certain 
QW j the field jumps down to either branches a or (3, and then decreases down towards 
the first fixed point E^^\r). Of course u has to be large enough for these jumps down 
to be possible. Thus these field profiles have a high-field domain located at the injecting 
contact and the domain boundary separating this domain from a low-field domain is a 
depletion layer. Numerical investigation shows that these stationary states are stable and 
that they can be reached from many initial conditions. Lemma 2 tells us that the initial 
field distribution can not be monotone increasing in the well index (like the solutions of the 
connected branch discussed before), for otherwise the field distribution would stay increasing 
for all times. Thus, these different stationary solutions are not connected to the branches 
dicussed before (having an increasing field profile) but form many additional isolated closed 
curves or "isolas" ||26| , |27[| . A typical isolated curve is shown inside the frame in Fig. |l] for 
i> = 1.0, which is also blown up to an enlarged scale for the sake of clarity. 

A special situation arises if c = holds. In this case the branches of increasing and 
decreasing field profiles become connected and there appears much additional degeneracy 
leading to an extremely complicated structure as observed in Ref. |jl5|. 



IV. STABILITY OF THE STATIONARY STATES 

Up to now we have only discussed the existence of stationary states, but not their sta- 
bility properties. First of all, several stability properties can be established by topological 
arguments PB| , |25| . If several branches overlap at a fixed voltage, each second branch has to 
be unstable by general reasons. For example the middle branch (exhibiting positive differ- 
ential conductivity) of the Z-shaped characteristics in Fig. |^(z/ = 0.15) has to be unstable. 
The remaining branches may exhibit further bifurcations. In order to elucidate this, we will 
perform a linear stability analysis for the states constructed in the previous section. By this 
method we will prove the following statements: 

A. For large doping No exceeding approximately e{Em — Em) / (q) ■ v~m/ (f m — "Vm) we find 
stable domain states. 

B. For very small products Nd{N — 1) the almost uniform states are also stable. 

For a medium range of doping in between these two limits we find self generated oscillations 



as reported in Ref. [10,11 1, and further discussed in a companion paper [23 



In order to perform the linear stability analysis, we substitute 



E,it) = E* + e^' (18) 
I{t) = r + e^'j (19) 



in (|^) and (^, thereby obtaining 

'v{E-) V 



Ae,^,-^:;X|£),,_^(e.-e..O. (20) 

eo = ei. (21) 



8 



These linear equations determine all as a function of A. The fixed bias condition 



Ee. = 0, (22) 

i=l 

then determines the possible eigenvalues A. For A 7^ we now introduce the variable 
Yi = Xci/j and the parameters 

b- - a. - ^ (23) 

^ - I , a, - ^ , [26) 

in Eqs. (|20|) and (^). Then we obtain: 

y _ ^ + or / A + ai + 6i = 1 ^24) 

* A + + 6i \ and Yi^i = -X/a., J ^ ' 

Within the scope of this linear stability analysis, we find that stationary states having 
all their fields Ei in the positive differential mobility region are stable, which coincides with 
our physical intuition. This is formulated in the following lemma: 

Lemma 3. If hi > Q holds for all i = 1, . . . N the real parts of all eigenvalues have to be 
negative, i.e., the state is stable. Furthermore the current-voltage characteristic exhibits a 
positive slope dl*/d^. 

We prove this Lemma by contradiction. Let us assume that Re{X) > holds with A 7^ 0. 
As 61 > we directly find that Re{Yi) > and — 1| < 1. In order to satisfy the voltage 
condition (|2^ ) we conclude that there must be at least one Yj with Re{Yj) < 0. This implies 
directly that also |Yj — 1| > 1 must hold. But we find: 

^, _ \ai(Yi_i - 1) -bi\ ^ \ai\\Yi_i - 1| + \bi\ . 

* \X + bi + ai\ \bi + ai\ 

Given that jli — 1| < 1, the last equation implies that — 1| < 1 for all i. Thus, the case 
Re{X) > 0, A 7^ is excluded. 
For A = we obtain 

biei = j and [h + ai)ei = j + atCi-i . (27) 



Therefore all Cj have the same phase as j and the voltage condition ( [24 ) cannot be satisfied 
unless Cj = 0, Vi and j = 0, which is the trivial case. In conlusion A = is not an eigenvalue. 
Furthermore A = describes the infinitesimal change along the curve of stationary states. 
Eq. (|27|) tells us, that j/ii > 0, Vi. Identifying dl* = j and dcj) = J2(^i we obtain a positive 
slope of the current- voltage characteristic, i.e., dl* /dip > 0.« 
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A. Stability for sufficiently large v 



Lemma 3 establishes that stationary field profiles with v'{E*) > 0, V? are linearly stable. 
These profiles include: (i) trivial ones where all the Ei belong to the same branch of v{E), 
and (ii) profiles where the negative differential mobility region is crossed in a single jump. 
This means that for a certain value E* < 1 [1 < j < N — 1) there exists -Eji+i > -Em 
with f{E*^^, I*) = E*. As f{E,I*) is increasing for E > E^, the necessary and sufficient 
condition for the existence of such a value E*_^-^ is f{Ejn, I*) < E*. This condition yields 

{Em - E*)v{Em) < y{P - Vm) (28) 

This inequality is first fulfilled for E* = 1 and /* = 1, which are the largest values of the 
respective quantities for the low field domain. This gives: 

. > (29) 

If this condition is fulfilled, domain states are possible which cross the negative differential 
mobility region in a single jump and must be stable therefore. 

Nevertheless there can be stable states even for smaller doping u, as Lemma 3 only yields 
a sufficient and not a necessary condition for stability. For our v{E) curve inequality (|29D 
yields u > 0.32. Indeed we do not find any self-sustained oscillations for u larger than the 
value u ~ 0.27, which is somewhat smaller than our estimation. Checking other v{E) curves 
we have always found oscillations up to a doping roughly 15 — 40% lower than determined 
by the bound (^). Thus, the bound is not only a sufficient condition but also a reliable 
rough estimate for doping above which the oscillations disappear. 

Transforming to physical units we obtain a surface charge density per well: 

No > , (30) 

(f M -Vm)q 

which should be a reasonable approximation for the necessary doping density. For the model 
equations from Prengel 0] we find N^, > 2.4 x 10^^/cm^. Actually, oscillations are found 



in the model up to roughly ^ 10^^/cm^ for the regular superlattice and for somewhat 



higher values for a slight amount of disorder [20 



If we regard domains with the high-field domain located at the injecting contact the 
same arguments yield a bound 

u > ^^^^ (31) 

for the existence of stable domains which is larger by a factor 1/vm- For our v{E) curve this 
corresponds to u > 0.97. Indeed we have found stable domains with an depletion layer for 
p = 1.0 as depicted in Fig. ^. This indicates that this type of domain only appears for larger 
doping. This might explain that two different locations of the high-field domain have been 
reported in the literature. In Ref. it is found to be located at the injecting contact for a 
superlattice with = 8.75 x 10^^ cm~^ while in Ref. ||2^ the high-field domain is located 
at the receiving contact for a different superlattice with = 1.5 x 10^^ cm~^. 
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B. Stability for sufficiently small v 



Now we want to show that for sufficiently small doping, the (connected) branches of 
stationary solutions are stable. Then no self-sustained oscillatory branches bifurcating from 
them can exist. In order to do this, we note that for very low voltages the stationary state is 
stable as indicated by Lemma 3: all field values of this state are in the range < -E* < 1. We 
now increase the voltage and study whether an instability may occur by checking whether 
it is possible to have \ = iuj with > for some 0. (The case A = yields the saddle-node 
bifurcation at the point with d(f)*/dl* = 0, which causes the switching to another branch of 
the Z-shaped characteristic but typically does not generate any oscillatory behavior.) In the 
appendix we show that this is possible only if 

rAT n ^ • {'^iVi-VCi) VlCi \ 

(i\ — Iju > mm < , — — — > (32) 

^ ^ 1 4ci ' 2C{I*-vi)c2j ^ ' 

holds, with 

vi = min \^Vm , (1 + c)'^] 

Cl = rma.XEi<E<Eh\d^nv{E) /dE\ 

C2 = rma.XE,<E<E,\d^lnv{E)/dE^\ 



C 



{N - l)i/ci 



/(7V-l).cA 

\ Vi- VCi I 



Here Ei,Eh denote the minimal and maximal values of the field for the stationary field 
profile. Note that for small u we find C ^ vi/{vi — uci) and furthermore the terms uci 
become negligible. Then the right side no longer depends on v and N but only on the shape 
of v{E) and the parameter c. 

If u is smaller and the inequality ( ^21) is violated, no bifurcating oscillatory branches can 
issue forth from the steady state which is thereby stable. 



The bound (|3^ is far too small due to the rough estimations made during its derivation. 
Therefore the number itself should not be used for quantitative investigations. Nevertheless 
we now have shown that the stationary states are stable for low doping and that in the limit 
of long superlattices the critical doping decreases as 



C. Consequences for the continuum limit 

With respect to the continuum limit, — oo, u ^ 0, L := Nv < oo we directly find 
that there exists a minimal length Lm such that the stationary state is linearly stable if 
L < Lm- This lower bound is given by the Eq. (|3^ with z/ = 0, C = 1 and it can be derived 
directly from the equations valid in the continuum limit, as we shall report elsewhere pO[ . 



Eq. (p2D and similar bounds derived for other boundary conditions constitute an explicit 



form of the well-known N^^L criterion of the Gunn effect The dimensionless length 

L (proportional to doping times the semiconductor length [^^) has to be larger than a 
certain number for the stationary solution to be unstable. 

Obviously the upper bound in the doping v (for the absence of the oscillatory regime) 
does not exist in the continuum limit (z/ — > 0, — > oo, L = Nv fixed). The discreteness is 
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essential for the field distribution to jump from the low-field region to the high field region 
without any fields exhibiting negative differential v{E) in between, (which stabilizes the 
field distribution). This explains that these stable stationary domains can not be found in 
the usual Gunn diode: for large enough dc voltage bias, the Gunn diode may have a stable 



stationary solution with a large field near the anode. p^j3l[] However inside the diode there 
will be a region where the field takes values on the branch of negative differential velocity, 
which is different from what happens in the SL. 

V. CONCLUSIONS 

In this paper we have shown how the complex stationary current-voltage characteristic 
exhibiting domain branches is generated continuously as the doping increases. For low 
doping the characteristic follows the local v{E) relation. If more charges are present, the 
characteristic becomes Z-shaped. When the doping is even larger, wiggles appear. For each 
doping the characteristic is connected, and the field profiles of all its different branches are 
monotone. The different disconnected branches observed experimentally correspond to the 
stable solution branches of the full stationary current-voltage characteristic. It would be 
very interesting to investigate whether it is possible to stabihze the unstable branches so 
that the full characteristic could be observed, as in the case of the double-barrier resonant- 
tunneling diode |35|j28[] . Additionally, for large doping there exist isolated branches on the 



full current-voltage characteristic having non-monotonic field profiles. 

The stability analysis shows that the almost uniform field profile is stable for low dop- 
ing. The critical doping above which time-periodic oscillations of the current appear is 
inversely proportional to the sample length for fixed superlattice parameters. This is the 
same situation as in the famous N^L criterion for the Gunn Diode. For yet larger doping the 
time-periodic oscillations of the current disappear: there is an upper critical doping above 
which there appear stable stationary solutions with two electric field domains (separated 
by an abrupt domain wall extending almost one period of the superlattice). Obviously, 
this is not possible for the conventional Gunn Diode due to the lack of discretization. It 
is important to mention that the upper critical doping needed to stabilize stationary do- 
main structures is higher for profiles having depletion layers instead of accumulation layers 
between the different domains. 

This paper gives (for the first time) analytical bounds of the interval of dimensionless 
doping (proportional to 2D doping over field at resonant tunneling peaks) on which self- 
sustained current oscillations may exist. To ascertain the influence of the SL parameters and 
the voltage bias on the current oscillations themselves, a detailed analysis of the dynamics of 
the model is necessary. This analysis together with experimental verification and predictions 
on the possibility of tuning the oscillation frequency with applied bias will be presented in 



a companion paper elsewhere. p3 
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APPENDIX A: PROOF OF STABILITY FOR SUFFICIENTLY SMALL v 

Here we prove that \ = iuj with c<j > can be an eigenvalue of the linearized system (pOD 
only if V is sufficiently large. 

In order to do this, we assume that a given stationary field profile {-Ef} exhibits A = 
iio with a; > and derive several necessary conditions for this. Restricting ourselves to 
increasing field profiles, must either be located on the first or second branch of the v{lT) 
curve, as otherwise the second branch is not reached which is a necessary condition for the 
instability according to Lemma 3. Let us now determine the smallest value Ei and the 
largest value E^ the stationary field profile \Ef\ may take. Ei is given by the value of E^ 
on the first branch for which the current takes on its minimal value, /* = 1, considering 
that the field must eventually take values on the second branch of v{E\ Eq. ([T5|) yields 
viEi) = 1/(1 + c) which determines the field Ei < 1. Eh is given by the largest value that 
E'^ can take on the third branch of v{E). Noticing that /* < 1 + c in Eq. (|15|), we can 
adopt Eh as the solution of v{Eh) = 1 + c from the third branch. Thus, Ei and Eh depend 
only on the v{E) curve and the parameter c but not on the field profile {-Ef}- For sake of 
convenience we introduce the following quantities: 

Vi -.= 111111 \vm,-^—\ (Al) 



1 + C 

Ci := /* maxEi<E<Eh 
C2 := /* maxEi<E<Eh 



d\nv{E) 



BE 
92 \nv{E) 



dE^ 



(A2) 
(A3) 



Then we have 



tti > — and < Ci Vi (A4) 



In the following we will assume that u is so small that the function f{E, I*) (/* fixed) is 
always invertible and furthermore 

Oj + bi > vi/u — ci > , Vi (A5) 

holds. 

For X = iuj, uj > 0, we have ReYi > 0, and |Yi — 1| < 1. As previously explained in the 
proof of Lemma 3, there must be a Yj such that KeYi > 0, i = 1, . . . , j — 1, and Rel^- < in 
order to fulfil the voltage condition. We are going to prove the following result: 
Lemma 4. Let j > 1 be the index that satisfies ReYi > 0, i = 1, . . . , j — 1, and ReYj < 0. 
(a) If u < ci, we have 
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AAci 

where A is the maximum of the expressions 

vi ^ 1 



0- - 1) - > 1^ (A6) 



A, 



Iw ~1. tti + bi 



{j -k + Vju ^ ai + h, 



fork = 2,.. .j. 

(b) If u > ci, we have 



where 



2B [I* - vi) C2 



B := , (A9) 



^^■■=max^.<E.<EA.^^^^^^. (AlO) 



Proo/ 

(a) Let oj < Ci\ 

In order to prove ([A6|) , we consider how the argument 0i of the complex quantity Yi is 
varying with i. 



Aeii4>^-.-4>.) = 1 (,^(1 _ Yr') + 6, + a^) (All) 



Therefore we get: 



Furthermore we have: 



i-i -'t>i = arctan — ■ ' ^ (A12) 

\hi + ai + ijjIm[Yi )) 



^1 = arctan I M (A13) 



Straightforward calculations starting from eq. (l24| ) yield: 

u? + ajCi;/m(Fj_i) + ai(6j + aj)-Re(Fj_i^ 



a;2 + (6i + 



2 



i?e(y,) = -^/m(y,) + ^^/2e(F,_i) (A14) 
Oj + fli Oi + ai 

By definition of the index j (j < A^), ReYj- < and ReY^_i > 0. Then these equations 
indicate that ImY, < and ImYj^i < 0. Thus the transition Rel^_i > ^ Rel^ < 
occurs across the angle = — 7r/2 as we have — 7r/2 < < and — vr < </)j < — 7r/2. 



14 



We introduce the index j' which is defined by the relations —n/A < (pj'-i and (pi < — 7r/4 
for i = j', . . . ,j. Obviously, we have ImFj < and therefore bi + ai — co^^f^ > for 
^ = + 1)---)J- Using ( A14 ) and R,eYi_i > (for all i < j), we obtain KeYi > uj 
\mYi/{ai + hi). This yields for i = j', . . . ,j: 



■>i-i 



arctan 



< arctan 



arctan 



Oi -\- tti — LlJ |y,|2 



UJ 



u)'^Im(Yi) 
(fei+a,)|y,|2 



Ol -\- Ul — UJ |y,|2 



a; 



< 



UJ 



tti + biJ ai + bi 



(A15) 



Now we have to distinguish two different cases: 

i) j' > 2: By summing the inequality ( |A15| ) from i = j' to i = j and then taking into 
account the definitions of j and /, we find: 



Vr TT ^ UJ 

uj{j -f + l)uAj, 



vi 



(A16) 



where definition (|A7|) has been used. The property uj < ci then implies 



J - 1 > J - j' + 1 > 



TTVi 



AAjivci 



(A17) 



ii) j' = 1: This means that 0i < —it/A and according to to Eq. ( |A13| ), uj < —hi = \bi\. Now 
we sum the inequality (|A15|) from i = 2to i = j and then use the expression arctan x > nx/A 
for < X < 1, thereby obtaining 



nuj 
-Ahi 



< arctan 



< 



UJ 



arctan 



ujv{j - 1) A2 



UJ 



< 



Vl 



TT 

2 



(A18) 



Therefore we find with Eq. 



J-l> 



TTVl 



> 



TTVl 



4Z/A2I&1I 4A2Z/C1 



Putting together (|A17|) and ( |A19| ) we obtain the inequality (|A6|) . 
(b) Let now u > Ci. 

]), we shall define the auxiliary functions 



To prove the inequality 



(A19) 



Zi — Yi 



X 

X + bi 



(A20) 
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These functions solve the following discrete equation 



(A + + bi) Zi - tti = Aflj 



\ + bi-i \ + biJ 

with the boundary condition Zi = 0. The solution of this problem is 

r7 _ A — bk-i) -A- tti 

" " h (A + (A + bk) M A + + a, • 



(A21) 



(A22) 



As all bi are real quantities, we have |zct; + 6j| > uj and obtain the following inequality for 
\Zj\ by using the preceding formula with A = iu, uj > 0: 



I j 3 

\Zj\ < - ^\bk- bk-i\ Y[ 



UJ 



fc=2 



bi + tti + icol 



(A23) 



Now we have \Zj\ > - ReZj = - ReYj + u;V(^^ + > t^V(^^^ + where ReYj < and 
the definition of Z„ have been used. This inequality together with ( A23|) yield 



UJ^ + b] 



3 3 
'k — bk~l 



k=2 



n — - — 

i=k h + «j + 



(A24) 



We now estimate the right side of ( |A24|) . The definition (^) of bi and the mean value 
theorem yield 



\bk - bk-i\ < C2 \E^k - Ek-i\ 



(A25) 



Equation ([T^) for the stationary state now yields < — E^_^ = {I*/ak — J^), so that 
< E^ — E^_^ < /vi — 1). Thus we can write: 



bk-i\ < V 



r 

Vl 



1 C2. 



On the other hand, as we are considering the case a; > ci, we find that 

di (Xi 



< 



\bi + ai + iuj\ \bi + ai + ici\ 
according to ( |A1CI| ). Inserting ( |A26| ) and ( [A27| ) into ( |A24| ) , we obtain 



(A26) 



(A27) 



UJ 



Uj^ + Cf 



Vl 



I/C2 I] fi 

Vl J yU - 1 



k=2 

/i^-l - 1 



(A28) 



Since uj > ci, uj^/{uj'^ + cf) > ci/2. Inserting this into (|A28|) , we obtain (|A8|) . Therefore 
Lemma 4 is proved.* 

Lemma 4 yields necessary conditions for the instability of a given stationary field profile 
{E'ji} corresponding to a fixed bias. We would like to obtain a general condition on z/, which 
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should only depend on the v{E) curve and the parameter c, but not on the specific stationary 
field profile. This can be achieved by the following considerations: 



At 



Vi - UCi 



(A29) 



Therefore we have A < vi/ {vi — vci ) from Eq. ( [A7[ ) , which inserted into the inequality (|A6| ) 
gives 



(j -l)y> IT 



4ci 



^From the definition (|A10D we obtain 

< meiXEi<E,<Eh 

This yields 



ai 



ai - ci 



< 1 + 



Z/Ci 



Vi - UCi 



B < 



< 



Vl 



1 + 



i-1 



exp 



to be inserted in 



). The result is 

(j -l)iy> 



Vl - Z/Ci . 

' {j - l)z/c/ 

, Vl - VCi , 



- 1 



VlCi 

2C {I* - Vl) C2 



(A30) 



(A31) 



(A32) 



(A33) 



We now use the obvious inequality > j in ( |A30D and ( |A33| ) thereby obtaining the 
condition (|32D as a necessary condition for oscillatory instability of the steady state. 



17 



REFERENCES 



* Now at: Mikroelektronik Centret, Technical University of Denmark, DK-2800 Lyngby, 
Denmark 

[1] L. Esaki and L. L. Chang, Phys. Rev. Lett. 33, 495 (1974). 

[2] Y. Kawamura, K. Wakita, H. Asahi, and K. Kurumada, Jpn. J. Appl. Phys 25, L928 
(1986). 

[3] K. K. Choi, B. F. Levine, R. J. Malik, J. Walker, and C. G. Bethea, Phys. Rev. B 35, 
4172 (1987). 

[4] P. Hclgcscn and T. G. Finstad, in Proc. l^th Nordic Semicond. Meeting (Arhus, 1990), 

ed. by Ole Hansen, p. 323. 
[5] H. T. Grahn, R. J. Haug, W. Miiller, and K. Ploog, Phys. Rev. Lett. 67, 1618 (1991). 
[6] S. H. Kwok, R. Merlin, H. T. Grahn, and K. Ploog, Phys. Rev. B 50, 2007 (1994). 
[7] J. Kastrup, H. T. Grahn, K. Ploog, F. Prengel, A. Wacker, and E. SchoU, 

Appl. Phys. Lett. 65, 1808 (1994). 
[8] R. Merhn, S. H. Kwok, T. B. Norris, H. T. Grahn, K. Ploog, L. L. Bonilla, J. Galan, J. A. 
Cuesta, F. C. Martinez, and J. M. Molera, in Proc. 22nd Int. Conf. Phys. Semicond., 
edited by D. J. Lockwood (World Scientific, Singapore, 1995), p. 1039. 
[9] S.-H. Kwok, T. C. Norris, L. L. Bonilla, J. Galan, J. A. Cuesta, F. C. Martinez, J. M. 
Molera, H. T. Grahn, K. Ploog and R. Merhn, Phys. Rev. B 51, 10171 (1995). 
[10] H.T. Grahn, J. Kastrup, K. Ploog, L. L. Bonilla, J. Galan, M. Kindelan, and M. 

Moscoso, Jpn. J. Appl. Phys. 34, 4526 (1995). 
[11] J. Kastrup, R. Klann, H.T. Grahn, K. Ploog, L.L. Bonilla, J. Galan, M. Kindelan, M. 

Moscoso, and R. Merhn, Phys. Rev. B 52, 13761 (1995). 
[12] A. A. Ignatov, V. I. Piskarev, and V. I. Shashkin, Sov. Phys. Semicond. 19, 1345 (1985). 
[13] B. Laikhtman, Phys. Rev. B 44, 11260 (1991). 

[14] F. Prengel, A. Wacker, and E. SchoU, Phys. Rev. B 50, 1705 (1994), erratum in 52, 
11518 (1995). 

[15] L. L. Bonilla, J. Galan, J. A. Cuesta, F. C. Martinez, and J. M. Molera, Phys. Rev. B 
50, 8644 (1994). 

[16] D. Miller and B. Laikhtman, Phys. Rev. B 50, 18426 (1994). 

[17] L. L. Bonilla, in Nonlinear Dynamics and Pattern Formation in Semiconductors, edited 

by F. J. Niedernostheide (Springer, Berlin, 1995), Chap. 1. 
[18] A. Wacker, F. Prengel, and E. SchoU, in Proc. 22nd Int. Conf. Phys. Semicond., edited 

by D. J. Lockwood (World Scientific, Singapore, 1995), p. 1075. 
[19] O. M. Bulashenko and L. L. Bonilla, Phys. Rev. B 52, 7849 (1995). 
[20] A. Wacker, G. Schwarz, F. Prengel, E. SchoU, J. Kastrup, and H.T. Grahn, Phys. Rev. B 

52, 13788 (1995). 

[21] A. N. Korotkov, D. V. Averin, and K. K. Likharev, Appl. Phys. Lett. 62, 3282 (1993). 
[22] A. N. Korotkov, D. V. Averin, and K. K. Likharev, Phys. Rev. B 49, 1915 (1994). 
[23] J. Kastrup, H. T. Grahn, R. Hey K. Ploog, L. L. Bonilla, M. Kindelan, M. Moscoso, 

A. Wacker and J. Galan, Phys. Rev. B 55, 2476 (1997). 
[24] O. M. Bulashenko, L. L. BoniUa, J. Galan, J. A. Cuesta, F. C. Martinez and J. M. 

Molera, in Quantum Transport in Ultrasmall Devices, edited by D.K. Ferry, H.L. Grubin, 

C. Jacoboni, and A.P. Jauho. NATO - AST Series B. (Plenum Press, New York, 1995), 

page 501. 



18 



M. Patra, Studienarbeit, Technische Universitat Berlin (1995). 

G. looss and D. D. Joseph, Elementary Stability and Bifurcation Theory. (Springer, New 
York, 1980). 

D. Dellwo, H. B. Keller, B. J. Matkowsky and E. L. Reiss, SIAM J. Appl. Math. 42, 
956 (1982). 

A. Wackcr and E. SchoU, J. Appl. Phys. 78, 7352 (1995). 

S. H. Kwok, U. Jahn, J. Menniger, H. T. Grahn, and K. Ploog, Appl. Phys. Lett. 66, 
2113 (1995). 

L. L. Bonilla, M. Kindelan, M. Moscoso and A. Wacker (unpubhshed) . 

M. P. Shaw, H. L. Grubin and P. R. Solomon, The Gunn-Hilsum Effect. (Academic 

Press, New York, 1979). 

F. J. Higuera and L. L. Bonilla, Physica D 57, 161 (1992). 

L. L. Bonilla, F. J. Higuera and S. Venakides, SIAM J. Appl. Math. 54, 1521 (1994). 
P. Gueret, Phys. Rev. Lett. 27, 256 (1971). 

A. D. Martin, M. L. F. Lerch, P. E. Simmonds, and L. Eaves, Appl. Phys. Lett. 64, 
1248 (1994). 



19 



FIGURES 



FIG. 1. Current- voltage characteristics for c = 0.01 and N = 20 and different values of v. The 
full line denotes the states where the electric field Ei is strictly monotone increasing in i. For 
ly = 1.0 there appear additional branches with non-monotonic field profiles Ei. They are isolated 
from the stationary branches having increasing field profiles. We have shown one such branch, 
which has also been blown up for the sake of clarity. The dotted line is the v{E) curve used 
throughout this paper. 

FIG. 2. f{E,I) and a trajectory Ei indicating decreasing i for various doping densities and 
currents, (a) = 0.15, / = 0.8. (b) z/ = 1.0, I = 0.8. (c) = 1.0, / = 1.008. (d) = 1.0, 
/ = 0.599 = Ic. 

FIG. 3. Different stationary electric field profiles for c = 0.01, v = 1.0, (p = 1.2, and N = 20. 
The crosses mark a state of the connected branch from Fig. l(i^ = 1.0) while the circles mark a 
state belonging to the isolated branch. 
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